Non-collinear magnetism of Cr nanostructures on Fe 3 ML/Cu(001): first— principles and 

experimental investigations 
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A combined experimental, using X-ray Magnetic Circular Dichroism and theoretical investigation, 
using full-potential Korringa-Kohn-Rostoker (KKR) Green function method, is carried out to study 
the spin structure of small magnetic Cr adatom-clusters on the surface of 3 monolayers of fee Fe 
deposited on Cu(OOl). The exchange interaction between the different Cr adatoms as well as between 
the Cr atoms and the Fe atoms is of antiferromagnetic nature and of comparable magnitude, leading 
due to frustration to complex non-collinear magnetic configurations. The presence of non-collinear 
magnetic coupling obtained by ab initio calculations is confirmed by the experimental results. 



The microscopic understanding of the magnetic prop- 
erties of small transition metal particles has recently at- 
tracted a lot of interest. Small clusters show magnetic 
behavior which is distinctly different from the respec- 
tive bulk material; for instance, even non-magnetic ma- 
terials can show magnetism as small clusters Q. Spin 
and orbital magnetic moments are usually significantly 
enhanced when the average coordination number of the 
atoms is reduced in a cluster Q. Of particular interest 
is the tendency of small clusters to favor non-collinear 
spin structures. In particular, for clusters on surfaces the 
formation of non-collinear spin structures is observed to 
relax possible spin frustration which results from a com- 
petition between intra-cluster and cluster-substrate inter- 
actions. A prototype of such a system are small Cr clus- 
ters on a ferromagnetic fee Fe/Cu(001) substrate. The 
latter one has already challenged experimentalists and 
theoreticians for several years and is generally acce pted 
to be ferromagnetic (FM) up 4 Fe monolayers |3l \4. Il2j. 
On the other hand Cr atoms tend to couple antiferro- 
magnetically to each other and lead in some cases 
to non-collinear magnetism. 

The aim of this work is to determine the complex mag- 
netism of Cr nanoclusters on Fe/Cu(001) using state-of- 
the-art theoretical and experimental techniques. Com- 
pared to bec Fe, where adatoms can approach each other 
at second-neighbor positions at the closest, an fcc(OOl) 
Fe substrate allows for a close-packed cluster geometry 
with first-neighbor positions between adatoms. This en- 
hances the intra-cluster exchange interactions. Further- 
more, compared to a Ni substrate, a Fe substrate has 
a much stronger exchange interaction with the adatoms 
and ad-clusters. These properties provoke magnetically 
frustrated states in antiferromagnetic ad-clusters (such 
as Cr) and make the particular combination of cluster 
atoms and substrate intriguing. 

We first present the ab-initio description of this sys- 
tem showing thus the importance of competing intra- 
cluster and cluster-substrate exchange interaction lead- 



TABLE I: Size and angles of magnetic moments in the calcu- 
lated clusters in the non-collinear configuration. The atoms 
separated by - have the same magnetic moments and rotation 
angles. 



Cluster Atom M (/x s ) 


<?n 


<M°) 


Trimer A, B-C 2.57, 2.92 


77, 156 


180, 


Tetra 1 A-D, B-C 2.50, 2.50 


111, 111 


0, 180 


Tetra 2 A, B-D, C 2.85, 2.87, 2.31 172, 176, 13 


0, 0, 180 


Penta A, B, C 2.44, 2.47, 2.17 
D, E 2.48, 2.89 


138, 85, 46 0, 180, 180 
155, 164 0, 



ing to non-collinear spin structures. The second step 
consists on depicting the measurements obtained with 
X-ray Magnetic Circular Dichroism (XMCD) in order to 
be compared with theory. 

— Theory. Our calculations are based on the Lo- 
cal Density Approximation (LDA) of density functional 
theory with the parametrization of Vosko et al.^ and 
use the full-potential non-collinear KKR-Green function 
methodpil Ha]. In the first step the surface Green func- 
tions are determined by the screened KKR method [Jjj 
for three monolayers of Fe (Fe3Mi) on Cu(OOl), using the 
LDA lattice parameter of Cu (~ 3.51 A). The obtained 
magnetic order of Fe is ferromagnetic, with moments for 
the surface layer of ms = 2.66/is and the two subsurface 
layers (mg_i = 2.09/j.b and ms-2 — 2.18/Zb). In the sec- 
ond step we calculate the Green function of the adatoms 
on the surface, by using the surface Green function as a 
reference. The resulting cluster of perturbed potentials 
includes the potentials of the Cr adatoms and the per- 
turbed potentials of several neighboring shells. Angular 
momenta up to / max = 3 were included in the expansion 
of the Green functions. We consider all impurities at the 
unrelaxed hollow position in the first vacuum layer. The 
orientations assigned to the spin moments of the impuri- 
ties are always relative to the orientation of the substrate 
moment, which we take as the global frame [l3j |. 

We start with single Cr adatoms, which are antiferro- 
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magnetically (AF) coupled to the substrate with an en- 
ergy difference AEfai-af = 564.9 meV to the ferromag- 
netic (FM) state, showing the strength of the interaction 
with Fc. The preference of the AF configuration shows up 
also in a considerably larger Cr local moment of 3.30/is 
compared to the metastable FM configuration (2.80/zg). 
Having established the strong AF coupling of the single 
adatom, we turn to the adatom dimers. We considered 
the two adatoms placed as first neighbors with differ- 
ent magnetic configurations. As it turns out the parallel 
orientation within the dimer, with both atoms antiferro- 
magnetically oriented to the substrate, is the most stable 
state with Cr moments of 3.02/is. Surprisingly, we did 
not find a non-collinear solution which would be expected 
by the following argument: Both Cr atoms tend to couple 
AF to the substrate as we have seen for single adatoms, 
but they also tend to couple AF to each other as expected 
for elements having a half filled e?-band and as explained 
by the Alexander- Anderson model0,^3|. These com- 
peting interactions, should lead to magnetic frustration 
and to a non-collinear solution. However, when we intro- 
duce a rotation of the moments in our ab-initio calcula- 
tions, both impurities relax back to the AF coupling to 
the substrate even if there is an AF between the dimer 
atoms. In order to explain this we introduce a Heisenberg 
model with the exchange parameters calculated from a fit 
to total energy calculations. We assume a classical spin 
Hamiltonian of the form H = — ^ Jij^iej, where e 

is a unit vector defining the direction of the magnetic mo- 
ment and i and j indicate the dimer atoms and their first 
Fe neighbors. Taking into account only first-neighbor in- 
teractions and neglecting the rotation of Fe moments, 
we rewrite the Hamiltonian for the dimer in the form: 

H = -Jcr-CrCOs(01 + 02)— 4 J Cr -Fe(cOS 01 + COS 2 )- The 

angle defining the non-collinear solution is obtained after 
a minimization of the Heisenberg Hamiltonian: cos(#i) = 

C0S(6> 2 ) = -2J Cr -Fc/</cr-Cr- If 2|J Cr -Fc| > |Jcr-Cr|, the 

angle is not defined and the non-collinear solution does 
not exist. This is realized in the present case: 2|Jo-Fc| = 
2 x 80.8meV > | J Cr -Cr| = 77.6 meV. Note that the Cr-Fe 
coupling constants are considerably smaller than for the 
single adatom. 

Following the same procedure as for the dimer we first 
investigated several collinear magnetic configurations for 
the most compact trimer on the surface, which has the 
shape of an isosceles rectangular triangle (see Fig. [2a)). 
After self-consistency all different collinear configurations 
converged to the FM solution with drastically reduced 
moments, i.e. 0.19/ab for atom A and 0.42/is for both 
atoms B and C. If we allow the directions of the moments 
to rotate, the moments strongly increase to "normal" val- 
ues (2.57 hb for atom A and 2.92/is for atoms B and C). 
As seen earlier, the AF interaction within the dimer was 
not high enough to compete with the adatom-substrate 
interaction. Addition of a third adatom to the system 
forces a rotation of the moments. The two second neigh- 





FIG. 1: (Color online) Complex magnetic states of Cr ad- 
clusters on Fe3ML/Cu(001) surface: (i) Side view of the trimer 
is shown in (a) with two Cr atoms pointing down (the sec- 
ond one cannot be seen), (ii) the front view of the T-shape 
tetramer (tetramer 2) is presented in (b), (iii) both sides of 
tetramer 1 are shown in (c) and (d), and finally the pentamer 
is shown in (e) and (f). 



boring impurities B and C have a moment tilted down 
by an angle of 156° and adatom A moment is tilted up 
with an angle of 77° (Table 1). Thus due to the re- 
duced Fe-Cr coupling the trimer atoms show a nearly 
antiferromagnetic configuration. The total energy differ- 
ences between the non-collinear and collinear solution is 
A-EncoI-fm = —360 meV/adatom. 

We extended our study to bigger clusters, namely 
tetramers and pentamer. Two types of tetramers were 
considered: tetramer 1 is the most compact and forms a 
square (Fig. ^d)), while tetramer 2 has a T-like shape 
(Fig. [IJb)). Concerning tetramer 1 (Fig. QJc) and (d)), 
two collinear configurations were obtained: the AF solu- 
tion, where all impurities are coupled antiferromagneti- 
cally to the substrate, and a Ferri solution where in each 
two neighbors within the tetramer are oriented AF to 
each other. The Ferri state is less stable than the AF 
as can be seen from the energy difference AEperri-AF = 
543.8 meV/adatom. The Ferri configuration is charac- 
terized by magnetic moments of —2.45 (1b and 2.29 /is 
whereas in the AF configuration the atoms carry lower 
moments (1.94 hb)- However, the ground state configu- 
ration is non-collinear (Fig.^c) and (d)) with an energy 
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difference AE'af-NcoI = 80 meV/adatom and a magnetic 
moment of 2.5 fiB carried by each impurity. One notices 
that the first neighboring adatoms are almost AF cou- 
pled to each other(the azimuthal angle <f> is either equal 
to 0° or to 180°) with all moments rotated by the angle 
9 = 111 . 

For tetramer 2 (see Fig. [D^b)) we obtained several 
collinear magnetic configurations. The most favorable 
one is characterized by an AF coupling of the three cor- 
ner atoms with the substrate. The moment of adatom 
C, surrounded by the remaining Cr impurities, is then 
forced to orient FM to the substrate. When we al- 
low for the direction of the magnetic moment to relax, 
we get a non-collinear solution having a similar picture, 
energetically close to the collinear one (A_E co i_Ncoi = 
2.3 meV/adatom). Adatom C has now a moment some- 
what tilted by 13° (/i = 2.31 /«b) whereas adatom A 
has a moment tilted in the opposite direction by 172° 
(fi = 2.85 /is). Adatoms B and D have a moment of 
2.87 hb with an angle of 176°. We note that tetramer 1 
with a higher number of first neighboring adatom bonds 
(four instead of three for tetramer 2) is the most sta- 
ble one (A£ ; te t2-tcti = 14.5 meV/adatom) with the non- 
collinear solution shown in Fig.QIc) and (d). 

To study the pentamer, we have chosen a structural 
configuration (Fig. me) an d (f)) with the highest number 
of first neighboring adatom bonds (five) . This pentamer 
consists on a tetramer of type 1 plus an adatom (E) and is 
characterized by a non-collinear ground state. Let us un- 
derstand the solution obtained in this case: Tetramer 1 is 
characterized by a non-collinear almost in-plane magnetic 
configuration (see Fig. QJe) and (f)). As we have seen, 
a single adatom is strongly AF coupled to the substrate. 
However, when one moves it closer to the tetramer it 
affects primarily the first neighboring impurity (adatom 
C) by tilting the magnetic moment from 111° to 46°. 
Adatom E is also affected by this perturbation and ex- 
periences a tilting of its moment from 180° to 164°. As a 
second effect, the second neighboring adatom, A, is also 
affected and suffers a moment rotation from 111 to 138°. 
The AF coupling between first neighboring adatoms is al- 
ways stable, thus adatom D has also a moment rotated 
opposite to the magnetization direction of the substrate 
with an angle of 155° (fx = 2.48 /is). As adatom B tends 
to couple AF to its neighboring Cr adatoms, its magnetic 
moment tilts into the positive direction with an angle of 
85°. 

— Experiment. The magnetic properties of the de- 
posited clusters have been determined experimentally 
using XMCD. The experiments were performed at the 
beamline UE56/1- PGM at the BESSY II storage ring 
in Berlin. The mass selected chromium clusters were 
generated using a UHV-cluster source JL5j and deposited 
in-situ onto ultrathin Fe layers epitaxially grown on a 
Cu(100) surface. The iron multilayers were prepared by 
evaporating iron from a high purity iron sheet onto the 
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FIG. 2: (Color online) Spin moment per number of d-holes 
per atom versus cluster size. Fig. a shows the experimental 
results for cluster-sizes up to 13 adatoms while Fig.b gives 
the comparison with ab initio results for cluster-sizes up to 
5 adatoms. Blue squares describe the experimental values, 
black diamonds show represent collinear configurations, red 
circles correspond to non-collinear configurations and a black 
line connects the ratios obtained in the magnetic ground 
states. 



clean copper crystal. Subsequently the iron films in the 
thickness range of ^3-5 monolayers (ML) were magne- 
tized perpendicular to the surface plane using a small 
coil. The magnetization of the iron films was monitored 
by recording Fe 2p XMCD spectra. Before cluster depo- 
sition Argon multilayers were frozen onto the iron surface 
at temperatures below 30K. A layer thickness of 10ML 
of argon was used to ensure soft landing conditions [lij . 
After depositing the clusters into the argon buffer lay- 
ers, the remaining argon was desorbed by flash heating 
the crystal to ~80K. Low sample temperatures in the 
range of 30K and low Cr coverages of 3% of a mono- 
layer (ML) were used to prevent cluster-cluster interac- 
tion. The measurements have been carried out at a base 
pressure p < 3 • 10~ 10 mbar. Spectra have been measured 
from the cluster samples prepared as described above, in 
a size range of 1 to 13 atoms per cluster. Every step of 
the preparation has been checked using X-ray Photoelec- 
tron Spectroscopy (XPS) and/or X-ray Absorption Spec- 
troscopy (XAS). The absorption signal has been mea- 
sured using the Total Electron Yield TEY, i.e. the sam- 
ple current. 

The experimental results for the magnetic moments 
per d-hole shown in Fig|5] have been determined apply- 
ing XMCD sum rules |l7J. After a standard background 
treatment, difference and sum spectra are generated by 
subtraction and addition of the x-ray absorption spectra 
measured with different photon helicity. A typical exam- 
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FIG. 3: (Color online) Absorption, difference and sum spectra 
of a Cr3 cluster. 



pie for a Cr3-cluster is shown in FigEJ For the application 
of the sum rules, the areas of the 2p 3 / 2 and 2pi/ 2 peaks 
of the difference spectra have been integrated. The 2p 3 / 2 
peak which shows positive and negative contributions in 
the difference spectra has been integrated over both con- 
tributions. The values for the spin and orbital magnetic 
moments (not shown) per d-hole have been corrected for 
a degree of circular polarization of 90%. The contribution 
of the spin magnetic dipole operator (T z ) in the values of 
the spin moments is presently ignored. 

We note that we have applied the XMCD sum rules 
without empirical corrections which have been proposed 
for chromium. We are well aware of the fact that due 
to the problems in the application of the XMCD sum 
rules for the early transition metals with smaller spin- 
orbit splitting the absolute values for the spin magnetic 
moments might be too small by up to a factor of two. 
However, this will not affect the relative trends seen 
in the experimental data (Fig. [3(a)). One notices the 
strong decrease of S e f f /rid with increasing cluster size 
which is due to the appearance of antiferromagnetic or 
non-collinear structures as calculated by theory. The 
qualitative and quantitative trends observed in the ex- 
perimental results agree very well with the theoretical 
results (Fig. EJ(b) for Cr-atoms to Cr-pentamers. Al- 
though the theoretical values for Cr-atoms and dimers lie 
somewhat higher than the experimental values (including 
error bars) the agreement can still be judged to be very 
good in view of the remaining experimental uncertainties 
discussed above. Even details as the increase of spin mag- 
netic moment from trimer to tetramer can be addressed. 
In order to understand the peak formed for the tetramer 
we compare in Fig. [2J (b) the ratio between the moment 
along the z-direction (defined by the magnetization of 
the substrate) and number of holes per atom obtained by 
theory for the different geometrical and spin structures. 
Black circles show the ratio calculated by taking into ac- 
count the collinear solutions and the red ones show the 



ratio calculated from the non-collinear solutions. The 
black line connects the ratio obtained in the magnetic 
ground states. The non-collinear tetramer 1 has a much 
lower value than what was seen experimentally whereas 
the collinear tetramer 1 and tetramer 2 give a better de- 
scription of the kink seen experimentally. With regards 
to the small energy difference (AE — 14.5V meV) be- 
tween the two tetramers we considered, we believe that 
the experimental value can be understood as resulting 
from an average of non-collinear tetramers 1, collinear 
and non-collinear tetramer 2. We believe that this ex- 
plains why the tetramer ratio value is higher than the 
one obtained for a trimer. The trimer and the pentamer 
are clearly well described by the theory and fit to the 
experimental measurements. 

To summarize, the competing exchange interactions 
in Cr-clusters and with Fe surface atoms affect strongly 
the complex magnetic structures leading thus to non- 
collinear magnetism in some cases. This is predicted 
by DFT calculations and confirmed by XMCD measure- 
ments. Considering the complexity of the different clus- 
ters and the non-collinear spin structures and the dif- 
ficult experimental analysis the quantitative agreement 
between theory and experiment is quite satisfying. A 
mechanism is proposed to explain the kink observed ex- 
perimentally for the tetramer in terms of geometrical con- 
figuration average. For further corroboration of this re- 
sult detailed STM investigations with spin-polarization 
analysis could be envisioned. 
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